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Abstract 

The use of the seniority quantum number in many-body systems is reviewed. A brief summary 
is given of its introduction by Racah in the context of atomic spectroscopy. Several extensions of 
Racah's original idea are discussed: seniority for identical nucleons in a single- j shell, its extension 
to the case of many, non-degenerate j shells and to systems with neutrons and protons. To illustrate 
its usefulness to this day, a recent application of seniority is presented in Bose-Einstein condensates 
of atoms with spin. 
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I. RACAH'S SENIORITY NUMBER 



The seniority quantum number was introduced by Racah for the classification of electrons 
in an atomic £ n configuration He assumed a spin- independent interaction V between the 
electrons with the property 

(£ 2 ;LM L \V\£ 2 ;LM L ) =g(2£+l)5 L0 , (1) 

that is, there is no interaction unless the two electrons' orbital angular momenta i are 
coupled to a combined angular momentum of L = 0. Racah was able to derive a closed 
formula for the interaction energy among n electrons and to prove that any eigenstate of 
the interaction ([1]) is characterized by a 'seniority number' v, a quantum number additional 
to the total orbital angular momentum L, the total spin 5* and the number of electrons n. 
He also showed that v corresponds to the number of electrons that are not in pairs coupled 



to L = 



a. 



Racah's original definition of seniority made use of coefficients of fractional 
parentage. He later noted that simplifications arose through the use of group theory jsj]. 
Seniority turned out to be a label associated with the orthogonal algebra SO(2£ + 1) in the 
classification 

U(4£ + 2) D (V(2£ + 1)D SO(2£ + 1) D ■ • O SO(3)) <g> SU S (2), (2) 

where the dots indicate intermediate algebras, if any exist. The number of states available to 
a single electron in an £ orbit is A£ + 2. All states of the £ n configuration therefore belong to 
the totally antisymmetric irreducible representation (IR) [l n ] of U(4£ + 2). Furthermore, the 
natural scheme for electrons in an atom is LS coupling which corresponds to the reduction 
U(4£ + 2) D U(2£ + 1) <g> SUs(2), where the orbital degrees of freedom are contained in 
U(2£+ 1) and the spin degrees of freedom in SUs(2). For any value of £ the unitary algebra 
XJ(2£ + 1) contains the orthogonal subalgebra SO(2£ + 1) which in turn contains SO(3), 
associated with the total orbital angular momentum L. 

The group-theoretical classification (J2]) allowed Racah to derive a number of important 
results in the theory of complex atomic spectra. The pairing force ([T]), however, is a poor 
approximation to the Coulomb interaction between electrons and for a more physically 
relevant application of seniority we have to turn to nuclei. 
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II. SENIORITY IN A SINGLE j SHELL 



The discussion of seniority in atoms and in nuclei differs in two aspects: (i) LS coupling 
is a good first-order approximation in atoms while in nuclei it is rather jj coupling and (ii) 
electrons are identical particles while nucleons come in two kinds, neutrons and protons. Let 
us postpone the discussion of the second complication until Sect. IIVI and concentrate in this 
section on the case of identical nucleons (either all neutrons or all protons). We impose the 
additional restriction that the identical nucleons are confined to a single- j shell, deferring 
the discussion of the many-j case to Sect. IIHI 

It turns out that a pairing force of the type 

(j 2 ; JMj\V\f; JMj) = -g(2j + 1)S J0 , (3) 

is a reasonable first-order approximation to the strong interaction between identical nucleons. 
In Eq. ([3]) j is the total (orbital+spin) angular momentum of a single nucleon and J results 
from the coupling of two of them. Since the pairing property now refers to the total j of the 
nucleons, there is no need for a separate treatment of orbital and spin degrees of freedom as 
in Eq. (T5]), and the classification becomes in fact simpler: 

U(2j + 1) D Sp(2j + 1) D ■■■ D SO(3) 

11 I • (4) 

[l n ] [l u ] J 

Seniority is associated with the (unitary) symplectic algebra Sp(2j + 1) which replaces the 
orthogonal algebra SO(2£+ 1) of the atomic case. Since the nucleons are identical, all states 
of the j n configuration belong to the totally antisymmetric IR [l n ] of U(2j + 1). The IRs of 
Sp(2j + 1) therefore must be totally antisymmetric of the type [1"]. The allowed values of 
seniority are v = n, n — 2, . . . , 1 or 0. The angular momentum content for a given seniority 
v can also be worked out Q] but no simple general rule is available. 

An alternative, simpler definition of seniority can be given which relies on the existence 
of an SU(2) symmetry of the pairing hamiltonian 5|, |6j. In second quantization the pairing 
interaction ([3]) is written as 

V = -gSiSL, (5) 

with 

S{ = ^27+1 (a] x a})?', Si = (stf , (6) 



where at creates a nucleoli in the shell j with projection mj. The commutator of S J + 
and Si. leads to the operator [S J + , SL] = (2fij — 2j — l)/2 = 2S J Z , which thus equals, up 
to a constant, the number operator hj. Since the three operators {S J Z ,S J ±} close under 
commutation, [S z , S 3 ±] = ±5*4 and [S+,S^.] = 2S Z , they form an SU(2) algebra, referred to 
as the quasi-spin algebra. 

This algebraic structure allows an analytical solution of the pairing hamiltonian. From 
the commutation relations it follows that S^SL = {Si) 2 — (Si) 2 + S£, which shows that 
the pairing hamiltonian can be written as a combination of Casimir operators belonging 
to SU(2) and SO(2) = {S J Z }. The associated eigenvalue problem can be solved instantly, 
yielding the energy expression — g[S(S + 1) — Ms(Ms — 1)]. The quantum numbers S and 
Ms can be put in relation to the seniority v and the nucleon number n, S — (2j — 2v + l)/4 
and Ms = (2n — 2j — l)/4, leading to the energy expression —gin — v)(2j — n — v + 3)/4. 
This coincides with the original expression given by Racah, Eq. (50) of Ref. after the 
replacement of the degeneracy in LS coupling, 4£ + 2, by the degeneracy in jj coupling, 
2j + 1. 

While this analysis shows that the eigenstates of a pairing interaction carry good seniority, 
it does not answer the question what are the necessary and sufficient conditions for a general 
interaction to conserve seniority. Let us specify a rotationally invariant two-body interaction 
V by the matrix elements vj = (j 2 ; JMj\V\j 2 \ JMj) with J — 0, 2, . . . , 2j — 1. The necessary 
and sufficient conditions for the conservation of seniority can then be written as 



^ , / I (3 3 4,/(2J+l)(2/+l) \ 

s V2JTT \s„ + 2V<2J + D(2/ + d{ ; ; f } - % _ +l) ' y^°> m 

with J = 2, 4, . . . , 2j — 1, and where the symbol between curly brackets is a Racah coefficient. 
These conditions have been derived previously in a variety of ways Q-Q. Although ([7]) 
determines all constraints on the matrix elements vj by varying / = 2, 4, . . . , 2j — 1, it does 
not tell us how many of those are independent. This number turns out to be [(2j — 3)/6j, 
the number of independent seniority v = 3 states [lOj. No condition on the matrix elements 
uj is obtained for j = 3/2, 5/2 and 7/2, one condition for j = 9/2, 11/2 and 13/2, and so 
on. As a result, identical nucleons in a single shell with j < 7/2 conserve seniority for any 
interaction 7]. 

Clearly, the conditions (JTj) are much weaker than the requirement that the interaction be 
of pairing character but still many of the results of the quasi-spin formalism remain valid. 



For instance, the ground state of an even-even nucleus still can be written in the form (JSJ). 
The main restriction of the concept of seniority as defined so far, concerns the fact that the 
nucleons are confined to a single- j shell. To lift this restriction, we turn to the generalization 
presented in the next section. 



III. SENIORITY IN SEVERAL j SHELLS 

The quasi-spin algebra can be generalized to the case of several degenerate shells (which 
we assume to be s in number) by making the substitutions S J + S + = J2j S+ and 2j + 1 i— > 
J2j{2j + 1)- Therefore, if a semi-magic nucleus can be approximated as a system of identical 
nucleons interacting through a pairing force and distributed over several degenerate shells, 
the formulas of the quasi-spin formalism should apply. In particular, the ground states of 
even-even semi-magic nuclei will have a 'superfluid' structure of the form 

n/2 



(s + no>, (8) 



where |o) represents the vacuum (i.e., the doubly-magic core nucleus). The SU(2) quasi-spin 
solution of the pairing hamiltonian (jSj) leads to several characteristic predictions: a constant 
excitation energy (independent of n) of the first-excited 2 + state in even-even isotopes, 
the linear variation of two-nucleon separation energies as a function of n, the odd-even 
staggering in nuclear binding energies, the enhancement of two-nucleon transfer. 

A more generally valid model is obtained if one imposes the following condition on the 
hamiltonian: 

[[H,S + },S + } = A(S + )\ (9) 

where S+ creates the lowest two-nucleon eigenstate of H and A is a constant. This condi- 
tion of generalized seniority, which was proposed by Talmi 11], is much weaker than the 
assumption of a pairing interaction and, in particular, it does not require the commutator 
[S + , S-] to yield (up to a constant) the number operator — a property which is central to 
the quasi-spin formalism. In spite of the absence of a closed algebraic structure, it is still 
possible to compute the exact ground-state eigenvalue but hamiltonians satisfying (jUj) are 
no longer necessarily completely solvable. 

An exact method to solve the problem of identical nucleons distributed over non- 
degenerate levels interacting through a pairing force was proposed a long time ago by 



Richardson [12J based on the Bethe ansatz [13[ . As an illustration of Richardson's approach, 



we supplement the pairing interaction with a one-body term, to obtain the following hamil- 
tonian: 

# = E eA - gs + s- = E -sE^iE^', ( 10 ) 

where 6j are single-particle energies. The solvability of the hamiltonian fflU]) arises as a 
result of the symmetry SU(2) ® SU(2) <g> • • • where each SU(2) algebra pertains to a specific 
j. Whether the solution of fflUj) can be called superfluid depends on the differences tj — 
6j> in relation to the strength g. In all cases the solution is known in closed form for 
all possible choices of ej. It is instructive to analyze first the case of n = 2 nucleons 
because it gives insight into the structure of the general problem. The two-nucleon, J = 
eigenstates can be written as S+\o) = J2j x jS+\o) with Xj coefficients that are determined 
from the eigenequation HS + \o) = ES + \o) where E is the unknown eigenenergy. With 
some elementary manipulations this can be converted into the secular equation 2ejXj — 
gJ2j> Qj>Xj' = Exj, with Qj = j + 1/2, from where Xj can be obtained up to a normalization 
constant, Xj oc g/(2ej — E). The eigenenergy E can be found by substituting the solution 
for Xj into the secular equation, leading to 

This equation can be solved graphically which is done in Fig. [1] for a particular choice of 
single-particle energies 6j and degeneracies Qj, appropriate for the tin isotopes with Z = 50 
protons and neutrons distributed over the 50-82 shell. In the limit g — > of weak pairing, the 
solutions E — > 2ej are obtained, as should be. Of more interest is the limit of strong pairing, 
g — > +oo. From the graphical solution we see that in this limit there is one eigenstate of the 
pairing hamiltonian which lies well below the other eigenstates with approximately constant 
amplitudes Xj since for that eigenstate \E\ 3> 2|ey|. Hence, in the limit of strong pairing one 
finds a J = ground state which can be approximated as 

^|o>«y^E^-l°>. ( 12 ) 

where Q = J2j Because of this property this state is often referred to as the collective S 
state, in the sense that all single-particle orbits contribute to its structure. 

This result can be generalized to n particles, albeit that the general solution is more 
complex. On the basis of the two-particle problem one may propose, for an even number 



y(E) 
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FIG. 1: Graphical solution of the Richardson equation for n = 2 fermions distributed over s = 
5 single-particle orbits. The sum ^ - Qj/(2ej — E) = y(-E') is plotted as a function of E\ the 
intersections of this curve with the line y = l/g (dots) then correspond to the solutions of the 
Richardson equation. 



of particles n, a ground state of the hamiltonian (ITUj) of the form (up to a normalization 
constant) 




E, 



Si |o), 



(13) 



which is known as the Bethe ansatz }13| . Each pair in the product is denned through 
coefficients Xj = (2ej — E a )~ x in terms of an energy E a depending on a which labels the n/2 
pairs. This product indeed turns out to be the ground state provided the E a are solutions 
of n/2 coupled, non-linear equations 



E 



E n 



n/2 

E 



Eft — E n 



a: 



l,...,n/2, 



(14) 



known as the Richardson equations 12] . Note the presence of a second term on the left- 
hand side with differences of the unknowns Ep — E a in the denominator, which is absent 
in the two-particle case. In addition, the energy of the state (ITBl is given by J2 a E a . A 
characteristic feature of the Bethe ansatz is that it no longer consists of a superposition of 
identical pairs since the coefficients (2e 3 - — E a )~ l vary as a runs from 1 to n/2. Richardson's 
model thus provides a solution that covers all possible hamiltonians ( !T0|) . ranging from those 



with superfluid character to those with little or no pairing correlations 

An important remaining restriction on the form of the pairing hamiltonian (JTU1) is that 
it contains a single strength parameter g whereas, in general, the interaction might depend 
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on j and j', leading to s(s + l)/2 strengths gjj> = g^j. In nuclei, often the assumption of a 
separable interaction is made which, in the case of pairing, leads to strengths gjj> = gCjCy 
in terms of s parameters Cj. This restriction leads to the following pairing hamiltonian: 

H = CjC r S 3 + S 3 l. (15) 

j jj' 

As yet, no closed solution of the general hamiltonian ffl5|) is known but three solvable cases 
have been worked out: 

1. The strengths Cj are constant (independent of j). This case was discussed above. 

2. The single-particle energies €j are constant (independent of j). The solution was given 



by Pan et al. [15] 



3. There are two levels. The solution was given by Balantekin and Pehlivan [161 ]. 

IV. SENIORITY WITH NEUTRONS AND PROTONS 

About ten years after its introduction by Racah, seniority was adop ted in nuclear physics 
for the jj-coupling classification of nucleons in a single- j shell [l8( . The main additional 
difficulty in nuclei is that one deals with a system of neutrons and protons, and hence 
the isospin T of the nucleons should be taken into account. The generalization of the 
classification (j3J) for identical nucleons toward neutrons and protons reads as follows: 

U(4j + 2) D (U(2j + 1) D Sp(2j + 1) D ■ ■ ■ D SO(3)) g> SU T (2) 

I I I 11, (16) 

[11 [h] [a] J T 

where [h] and [a] are Young tableaux associated with U(2j + 1) and Sp(2j + 1). In general, 
2j + 1 labels are needed to characterize an IR of U(2j + 1), [h] = [hi, h 2 , . . . , h 2 j + i], and 
j + 1/2 labels are needed for an IR of Sp(2j + 1), [a] = [a%, a 2 , ■ ■ ■ , crj + i/ 2 ]. To ensure 
overall antisymmetry under U(4j + 2), the Young tableaux of U(2j + 1) and Ur(2) must be 
conjugate, that is, one is obtained from the other by interchanging rows and columns. Since 
the Young tableau associated with Ux(2) is determined by the nucleon number n and the 
total isospin T as [n/2 + T, n/2 — T], the Young tableau of U(2j + 1) must therefore be 

n/2-T 2T 

[h] = f2,2,r..,2,l,l,...,l]. (17) 



Since an IR of U(2j + 1) has at most 2j + 1 labels, it follows that n/2 + T < 2j + 1. 
Furthermore, all non-zero labels in [a] must be either 2 or 1 and the Young tableau of 
Sp(2j + 1) must therefore be of the form 

v/2-t 2t 

[a] = f2,2,r.. s 2,i,l,r ..,f). (18) 

The IR of Sp(2j + 1) is thus characterized by two labels: the seniority v and the 'reduced 
isospin' t. The former has the same interpretation as in the like-nucleon case while the latter 
corresponds to the isospin of the nucleons which are not in pairs coupled to J = 0. 

The group-theoretical analysis is considerably more complex here than in the case of 
identical nucleons and, in addition, for each value of j one is faced with a different reduction 
problem associated with U(2j + 1) D Sp(2j + 1) D SO(3). It is therefore advantageous to 



go over to a quasi-spin formulation of the problem and, as was shown by Helmers |6|, this 
is possible for whatever value of the intrinsic quantum number of the particles (which is 
t = 1/2 for nucleons). If the pairing interaction is assumed to be isospin invariant, it is the 
same in the three T = 1 channels, neutron-neutron, neutron-proton and proton-proton, 
and Eq. ([5]) can be generalized to 

V' = -gJ2S+J-,» = -gS + ■ (19) 
n 

where the dot indicates a scalar product in isospin. In terms of the nucleon creation operators 
a jm tmv wri ich now carry also isospin indices (with t = 1/2), the pair operators are 

S + ,» = W \f^+~Mt x 4C> = (§ + jf , (20) 

L 3 

where the coupling refers to angular momentum and to isospin. The index /i (isospin projec- 
tion) distinguishes neutron-neutron (/i = +1), neutron-proton (/i = 0) and proton-proton 
(fi = — 1) pairs. There are thus three different pairs with J = and T = 1 and they are 
related through the action of the isospin raising and lowering operators T± . By considering 
the commutation relations between the different operators, a closed algebraic structure is 
obtained, generated by the pair operators S± t/i , the number operator n and the isospin oper- 



ators T± and T z . The quasi-spin algebra of neutrons and protons in degenerate j shel 
out to be SO (5), by virtue of which the hamiltonian f ll9p is analytically solvable 19 



s turns 



20|. 



A further generalization is possible in LS coupling. For a neutron and a proton there 
exists a different paired state with parallel spins. The most general pairing interaction for a 



system of neutrons and protons is therefore of the form 



V" = -gS + 'S.-g'P + -P-, 



(21) 



where the pair operators are defined as 



= ]fl E v / 27TT(aJ 



ist x a lst ) OOAt 5 



x a 1 " ) (010) 
est x a ^st ) O/Lto ) 



(22) 



where a 



im e sm s tm t crea t es a nucleon in the shell £ with projection m £ , spin projection m s and 
isospin projection m t . The hamiltonian (12T|) contains two parameters g and g', the strengths 
of the isovector and isoscalar components of the pairing interaction. While in the previous 
case the single strength parameter g just defines an overall scale, this is no longer true for a 
generalized pairing interaction and different solutions are obtained for different ratios g/g'. 
In general, the eigenproblem associated with the interaction 02 II) can only be solved 



numerically; for specific choices of g and g' the solution of V" can be obtained analyti- 



cally 



21 



n 



|22| . A closed algebraic structure is obtained, formed by the pair operators (1221) . 
their commutators, the commutators of these among themselves, and so on until closure is 
attained. The quasi-spin algebra in this case turns out to be SO (8), with 28 generators, 
consisting of the pair operators S± t/i and -P±, M , the number operator h, the spin and isospin 
operators and T M , and the Gamow-Teller-like operator Y^ u , which is a vector in spin and 
isospin. The symmetry character of the hamiltonian ( 121]) is obtained by studying the subal- 
gebras of SO(8). Of relevance are the subalgebras SOt(5) = {S±^,h,T^}, SO^(3) = {T M }, 
S0 5 (5) = {P± )M , n, S^}, SO s (3) = {S^} and SO(6) = {S^f^Y^}, which can be placed in 
the following lattice of algebras: 



> D S0 5 (3)®SO T (3). 



(23) 



SO s (5) g> SO T (3) 
SO(8) D I SO(6) 

^ SO T (5)®SO s (3) ( 

By use of the explicit form of the generators of SO (8) and its subalgebras, and their com- 



mutation relations 



22J, the following relations can be shown to hold: 

1 



S + -S_ = ^(7 2 [SO T (5)] - i(7 2 [SO T (3)] 



(2fi-n)(2fi-n + 6), 
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S + ■ + P + ■ P_ = ic 2 [SO(8)] - JOs[SO(6)] - ^(20 - n)(2fi - n + 12), 

2 2 8 

P+-P- = \c 2 [$O s {h)\ - i<%[SO s (3)] - i(2fi - n)(2fi - n + 6), (24) 

with f2 = J2g(2£ + 1) and where C n [G] is the n th -order Casimir operator of the algebra G. 
This shows that the interaction (121 j) in the three cases (i) g = 0, (ii) (?' = and (iii) g = (7', 
can be written as a combination of Casimir operators of algebras belonging to a chain of 
nested algebras of the lattice (T23]) . They are thus the dynamical symmetries of the SO (8) 
model. 

The nature of £ SO(8) superfluidity' can be illustrated in the specific example of the 
ground state of even-even N = Z nuclei. In the SO (6) limit of the SO (8) model the exact 
ground-state solution can be written as 

(S + -S + -P + -P + ) n/A \o). (25) 

This shows that the superfluid solution acquires a quartet structure in the sense that it 
reduces to a condensate of bosons each of which corresponds to four nucleons. Since the 
boson in f!25jl is a scalar in spin and isospin, it can be thought of as an a particle; its orbital 
character, however, might be different from that of an actual a particle. A quartet structure 
is also present in the two SO (5) limits of the SO (8) model, which yields a ground-state wave 
function of the type (1231) with either the first or the second term suppressed. A reasonable 
ansatz for the N = Z ground-state wave function of the SO (8) pairing interaction (|2"TT) with 
arbitrary strengths g and g' is therefore 

(costf S + ■ S + - sin6 P + ■ P + ) n/A |o), (26) 

where 6 is a parameter that depends on the ratio g/g'. The condensate ( 1261) of a- like 
particles provides an excellent approximation to the N = Z ground state of the pairing 
hamiltonian (l2Tj) for any combination of g and g' [3]. It should nevertheless be stressed 
that, in the presence of both neutrons and protons in the valence shell, the pairing hamilto- 
nian (1211) is not a good approximation to a realistic shell-model hamiltonian which contains 
an important quadrupole component. 

These results can be generalized to the case of several non-degenerate shells. In fact, the 



Richardson equations 



are valid for the quasi-spin symmetry SU(2) but they are known 



for any Lie algebra [241 ] . Closed solutions have been obtained for a system of neutron and 

11 



protons with a pairing interaction of pure isovector character and of equal isovector and 
isoscalar strength, based on the SO (5) and the SO (6) quasi-spin algebras, respectively [25|, 
26|- 



V. BOSE EINSTEIN CONDENSATES OF ATOMS WITH SPIN 

In this section the concept of seniority is illustrated with an application to the physics of 
cold atoms. If atoms in a Bose-Einstein condensate (BEC) are trapped by optical means j^ . 
their hyperfme spins (or spins) are not frozen in one particular direction but are essentially 
free but for their mutual interactions. As a result, the atoms do not behave as scalar particles 
but each of the components of the spin is involved in the formation of the BEC. This raises 
interesting questions concerning the structure of the condensate and how it depends on the 
spin-exchange interactions between the atoms. 



Such questions were addressed in a series of theoretical papers by Ho and co-workers [28 ] 
who obtained solutions based on a generating function method. In the case of spin-1 atoms 
the problem of quantum spin mixing was analyzed by Law et al. [29) who proposed an 
elegant solution based on algebraic methods. It is shown here that an exact solution is 
also available for the spin value / = 2 (for any number of atoms n) which allows the 
analytic determination of the structure of the ground state of the condensate. This was 



simultaneously and independently pointed out in Refs. [30|, |31|. 

We consider a one-component dilute gas of trapped bosonic atoms with arbitrary (integer) 
hyperfme spin /. In second quantization the hamiltonian of this system has a one-body and 
a two-body piece that can be written as (in the notation of Ref. [29]) 

where h — 1, M a is the mass of the atom, and i£? m and are the atomic field anni- 
hilation and creation operators associated with atoms in the hyperfme state \fm) with 
m = —f, . . . , +/, the possible values of all summation indices in (12"7|) . The trapping poten- 
tial Vtrap is assumed to be the same for all 2f + 1 components. According to the assumptions 
outlined in Ref. j^J, the atomic field creation and annihilation operators at zero tempera- 
ture can be approximated by W m w b^cp^x), ifr m ~ b m <p(x), m = —f, . . . , +/, where <f)(x) 
is a single wave function (independent of m) and b m and b^ are annihilation and creation 

12 



operators, satisfying the usual boson commutation rules. In this approximation the entire 
hamiltonian f )27|) can be rewritten as 

V, ^ H = eb ] -b+-J2 x 6 f ) (F) ■ (b x b) {F \ (28) 

2 F 

where the coefficients e and vp are related to those in the original hamiltonian ( 12?]) and with 

u m — V / u —m- 

Exactly solvable hamiltonians with rotational or SO (3) invariance are now found by the 
determination of all Lie algebras G satisfying U(2/ + 1) D G D SO(3). The canonical 
reduction of U(2/ + 1) is of the form as encountered by Racah (see Sect. 

U(2/ + l) D SO(2/ + l) D SO(3), (29) 

defining a class of solvable hamiltonians of the type 

H> = a x C x [\J{2f + 1)] + a 2 C 2 [XJ(2f + 1)] + bC 2 [SO(2f + 1)] + cC 2 [SO(3)], (30) 

where a%, a 2 , b, and c are numerical coefficients. The solvability properties of the original 
hamiltonian ( 128]) now follow from a simple counting argument. For atoms with spin / = 1 
the solvable hamiltonian ( 1301) has three coefficients a\, a 2 , and c [since SO(2/ + l)=SO(3)] 
while the general hamiltonian (128]) also contains three coefficients e, Uq, and v 2 . (Note that 
the coupling of two spins to odd F is not allowed in the approximation of a common spatial 
wave function, so no v\ term occurs.) They can be put into one-to-one correspondence. 
For atoms with spin f — 2 both the solvable and the general hamiltonian contain four 
coefficients (ai, a 2 , b, and c versus e, h> , u 2 , and z/ 4 ) which also can be put into one-to- 
one correspondence. Hence the general hamiltonian (j28~I) is solvable for / < 2. The same 
counting argument shows that it is no longer solvable for / > 2. 

The case of interacting / = 1 atoms was discussed by Law et al. [29] who identified the 
existence of two possible condensate ground states: one with all atoms aligned to maximum 
spin F = n and a second with pairs of atoms coupled to F = 0. Whether the condensate is 
aligned or paired depends on a single interaction parameter. With the technique explained 
above, the phase diagram for atoms with spin f — 2 can also be derived. The results 
are exact and valid for arbitrary n. The entire spectrum is determined by the eigenvalue 
expression together with the necessary branching rules. In particular, the allowed values of 
total spin F for a given seniority v are derived from the SO (5) D SO (3) branching rule given 
by F — 2t, 2t — 2, 2r — 3, . . . , r + 1, r with r = v, v — 3, v — 6, . . . and r > 0. 
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It is now possible to determine all possib 
and their quantum numbers Vq and Fq 



e ground-state configurations of the condensate 
30|. The character of the ground state does not 
depend on the coefficients a, t since the first two terms in the expression ( I30p give a con- 
stant contribution to the energy of all states. Although this contribution is dominant, the 
spectrum-generating perturbation of the hamiltonian is confined to the last two terms and 
depends solely on the coefficients b and c which are related to the original interactions vp 
according to b = (—7u + KV2 — 3^) /70 and c = (—^2 + v^j/lA. The phase diagram displays 
a richer structure than in the / = 1 case. There is an aligned phase where the seniority 
is maximal, Vq = n, and all spins are aligned, F = 2n. Secondly, there is a low-seniority 
(paired) and consequently low-spin phase. For even n, this corresponds to (vq,F ) = (0,0). 
The aligned and paired phases are also encountered for interacting / = 1 atoms. For f — 2 
a third phase occurs characterized by high seniority (i.e., unpaired) and low total spin, 
(i; 0> F Q ) = (n, 25) with 5 = or 1. 

Since the hamiltonian ( |30|) is solvable for / = 2, all eigenstates, and in particular the 
three different ground states, can be determined analytically. The general expressions given 



\v = 
\v 



by Chacon et al. [32| reduce to 

n,F — Mp = 2n) oc (^ 2 )"|0), 
= 0,F = M F = 0) oc (rf t -d t ) n/2 |0), 
\ v = n , F = M F = 0) oc (V x a f ) (2) • a f ) n/3 |0), (31) 

where the / = 2 atoms are denoted as d bosons. In the second of these expressions it is 
assumed that n is even and in the third that n = 3k; other cases are obtained by adding a 



single boson d) or a d) ■ d) pair. The a) are the so-called traceless boson operators 32] which 



are defined as = d^ — d) ■ d^(2n + 5) _1 c? m (see also Chapt. 8 of Ref. |33[]). The wave 
functions fl3Tl) are the exact finite-n expressions for the eigenstates of the hamiltonian (1301) . 
Since in the large-n limit the traceless boson operators become identical to d^, one 
arrives at a simple interpretation of the three types of configurations: (i) spin-aligned, (ii) 
condensed into pairs of atoms coupled to F = 0, and (iii) condensed into triplets of atoms 
coupled to F — 0. 

In conclusion, the consideration of seniority is crucial in obtaining results concerning 
Bose-Einstein condensates consisting of atoms with spin. Since all eigenstates of interacting 
atoms with spin / < 2 are known analytically, this opens up the possibility to study the 
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relaxation properties of such condensates using their exact, macroscopic wave functions. In 
addition, preliminary studies indicate that seniority can be exploited even when / > 2. 
These problems are currently under investigation (341 ] . 

This paper is dedicated to the memory of Marcos Moshinsky. The two years I have spent 
in Mexico as a visitor and the many hours with Marcos as a teacher, were crucial to my 
formation as a physicist. Without him I never could have written this paper. 
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